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We propose a phenomenological study of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) approach 
applied to the data on the proton structure function F2 measured at HERA in the small- x Bj 
region. In a first part we use a simplified "effective kernel" approximation leading to few-parameter 
fits of F2. It allows for a comparison between leading- logs (LO) and next-to-leading logs (NLO) 
BFKL approaches in the saddle-point approximation, using known resummed NLO-BFKL kernels. 
The NLO fits give a qualitatively satisfactory account of the running coupling constant effect but 
jy-^ ] quantitatively the \ 2 remains sizeably higher than the LO fit at fixed coupling. In a second part, 

, a comparison of theory and data through a detailed analysis in Mellin space xsj —> w, leads to a 

■ more model independent approach to the resummed NLO-BFKL kernels we consider and points out 

' some necessary improvements of the extrapolation at higher orders. 

a 

! I. INTRODUCTION 

. Precise phenomenological tests of QCD evolution equations are one of the main goals of deep inelastic scattering 
' phenomenology. For the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution in Q 2 [1], it has been possible 
to test it in various ways with NLO (next-to- leading logQ 2 ) and now NNLO accuracy and it works quite well in a large 
CO ■ range of Q 2 and XBj- Testing precisely the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution in energy [2] (or xsj) 
beyond leading order appears more difficult. A theoretically convenient way would be to stay within the perturbative 
regime by using only massive or highly virtual colliding particles, but precision QCD phenomenology for the present 
I day is provided mainly by the data [3] on deep-inelastic scattering at small xb j ■ 

Indeed, the first experimental results from HERA confirmed the existence of a strong rise of the proton structure 
function F2 with energy which, in the BFKL framework, can be well described by a simple (3 parameters) LO-BFKL 
fit [4, 5] . The main issue of Ref. [4] was that not only the rise with energy but also the scaling violations observed 
at small xsj are encoded in the BFKL framework through the Q 2 variation of the effective anomalous dimension. 
Oh, However one was led [4] to introduce an effective but unphysical value of the strong coupling constant 

a ~ .07— .09 instead of a ~ .2 in the Q 2 -range considered for HERA small- x Bj physics, revealing the need for NLO 
corrections. Indeed, the running of the strong coupling constant is not taken into account. 

In fact, the theoretical task of computing these corrections appears to be quite hard. It is now in good progress but 
still under completion. For the BFKL kernel, they have been calculated after much efforts [6]. In fact, they appeared 
to be so large that they miss by a large amount the phenomelogical requirements and could even invalidate the whole 
theoretical approach. Soon after, it was realized [7] that the main problem comes from the existence of spurious 
singularities brought together with the NLO corrections, which ought to be cancelled by an appropriate resummation 
at all orders of the perturbative expansion, resummation required by consistency with the QCD renormalization 
group. 

Indeed, various resummation schemes have been proposed [7-9] which satisfy the renormalization group require- 
ments while retaining the computed value of the NLO terms in the BFKL kernel. Hence, the constraints can be 
satisfied and the next-to-leading order introduced without compromising the theoretical consistency of the BFKL 
scheme. However, the situation remains not so clear concerning phenomenology 1 . 

To summarize the theoretical problems still remaining to be solved, the determination of the impact factors associ- 
ated with the coupling of the NLO-BFKL kernels with the virtual photon are still in progress [12]. The factorization 
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1 There exists fruitful phenomenological approaches starting from the DGLAP evolution and adding log xgj correction terms in the 
perturbative expansion [10, 11]. 
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of the non perturbative coupling is also problematic at NLO-BFKL level. Another source of indeterminacy comes 
from instabilities at large energy of the evolution equations which may lead to the dominance of the non-perturbative 
"Pomeron" singularity [13]. 

On a phenomenological ground, we note that the rcsummation schemes possess some ambiguity, since higher 
order logs (beyond the next-to-leading ones) are not known, apart from the renormalization group constraints. Such 
variations appear, e.g. in Ref. [7], where first resummation schemes 2 have been proposed. In practice, we will consider 
the schemes S3 and 54 of Ref. [7] together with the one of Ref. [8], in the formulation of Ref. [14], including quark 
contributions to the kernel. 

The aim of our paper is to compare the phenomenological efficiency of LO versus NLO BFKL approaches using 
present precise experimental data on the proton structure function F 2 . On the same footing it is also to compare 
phenomenologically the different NLO schemes in order to check their validity and distinguish between different 
resummation options for the NLO-BFKL kernels. 

For this sake, we shall use a two-step approach. First, we will consider a simpler version of the NLO-BFKL formulae 
by considering an "effective kernel" obtained using a well-known "consistency condition" [8] satisfied by the NLO- 
BFKL kernels. Then, we formulate a saddle-point approximation for both cases allowing a similar formulation for LO 
and NLO kernels 3 . 

Our study is organized as follows. In Sec. II, we present the construction of the effective kernels at LO and NLO 
levels and the corresponding saddle point approximation for the proton structure function F 2 . In Sec. Ill, we use these 
formulations to perform fits to the structure function data in the range (x B j < .01 , Q 2 < 150 GeV 2 ) suitable for a 
BFKL analysis. In Section IV we perform a phenomenological determination of the Mcllin transform F2(u>, Q 2 ) in the 
relevant range and, with this information, we discuss the "consistency condition" for the NLO-kernels under study. 
The final Section V is devoted to a summary, a discussion of our approximations and an outlook. An Appendix recalls 
the necessary formulae and definitions of the NLO kernels considered in this work. 



II. "EFFECTIVE KERNEL" AND SADDLE-POINT APPROXIMATION OF BFKL AMPLITUDES 

The BFKL formulation of the proton structure functions can be formulated in terms of the double inverse Mcllin 
integral 



^=//!$(S)V^>. 



At LO level one has (see, e.g. [4]) 



w - a XLoh) 

where a = a s N c /-K, a s is the coupling constant which is merely a parameter at this LO level. The LO BFKL kernel 
is written as 

X lo(7) = 2V(1)-V'(7)-V'(1-7) • (3) 

/i2( 7 ,c<j) is a prcfactor which takes into account both the phenomenological non-perturbative coupling to the proton 
and the perturbative coupling to the virtual photon. Note that the variable 7 plays the role of a continuous anomalous 
dimension while to is the continuous index of the Mellin moment conjugated with the rapidity Y = log 1/xsj- 

Recalling well-known properties of LO-BFKL amplitudes, one assumes that /i2(7,w) is regular. The pole contribu- 
tion at oj — a xwij) in (2) leads to a single Mellin transform in 7 for which one may use a saddle-point approximation 
at small values of XBj- Indeed, in the LO problem, it is known that the saddle-point approximation gives a very good 
account of the phenomenology as we will confirm later on. Going beyond the saddle-point approximation is theoreti- 
cally more accurate, but leads to the same quality of fits, see, e.g. [15]. 

The saddle-point approximation gives 

F ^ Q2) ™ N CXP { \ + - 2 a f xU l) } ■ (4) 



2 Other schemes have been proposed, e.g. [9], which are left for further study. 

3 A more complete analysis will be possible within the same framework when the full NLO BFKL analysis including NLO impact factors 
will be available . 
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where L = log(Q 2 /Qq) and M is a normalisation taking into account all the smooth prefactors 4 . As a consequence, 
the only three relevant parameters in (4) are TV, a and Qq. In this picture a has to be considered as a parameter 
and not a genuine QCD coupling constant since the value obtained in the fits is not related to the coupling constant 
values in the considered range of Q 2 . As we shall now see, an effective saddle-point expression similar to (4) for the 
NLO BFKL analysis can be written, but it will retain the running property of the QCD coupling constant with its 
theoretically predetermined value at the relevant Q 2 range. Hence it is no more a free parameter. 

We will assume that, for high enough virtuality Q 2 , NLO-BFKL solution for small-x^j structure function is dom- 
inated by the perturbative Green function. For this Green function [7, 8], a consistency condition relation holds, 
namely 

XNLo{j,u) 2 

w — = ui - a RG (Q ) xnlo{i, w) = , (5) 

where 

[anciQ 2 )]' 1 log (Q 2 / A 2 QCD ) , (6) 

with b = 11/12 - 1/6 N f /N c . 

The implicit relation (5) can be considered in two different ways. On the one hand, considering it as an implicit 
equation keeping ui,Q 2 fixed defining 7(w, a RG ), it corresponds to the saddle-point in the integration of the Green 
function over 7. On the other hand, keeping 7, Q 2 fixed, it can be considered as an implicit equation for a; (7, a RG ), 
and appears as the NLO counterpart of the pole dominance from (2). Starting with the relation (5), one defines 5 an 
effective NLO BFKL kernel 

Xeff(l,a RG ) = • (7) 

CtRG 

Using this kernel, and performing a saddle-point approximation at large Y on the amplitude, we obtain 

r2 



F 2 (x, Q 2 ) « N exp < 7c L + a RG Xeffilc, a RG ) Y - — j — 

{ 2a RG x" f f (lc, org) Y 



(8) 



where j c is defined by the implicit saddle-point equation 

d Xeff 



c?7 



(7c«fl G (Q 2 )) = 0. (9) 



It is important at this stage to notice that the formula (8) has only two free parameters M and Qq instead of three for 
(4), once using the QCD universal coupling constant a RG . It allows one to compare in a similar footing the LO and 
NLO BFKL kernels to F 2 data. In first place, it allows for examining the important effect of the running coupling 
constant on the fits. 



III. F 2 FITS USING LO AND NLO BFKL KERNELS 



In practice, the determination of the effective kernel Xeffil, &rg) using the implicit relation (7) is made with the 
input corresponding to the (rcsummcd) NLO BFKL kernels proposed in the literature. As recalled in the introduction, 
the NLO BFKL kernels have to be properly defined, in order to incorporate the next-leading terms calculated in [6] 
and to get read of spurious singularities which would contradict the renormalization constraints [7-9]. There are 
different options for satisfying these constraints. As an input of our analysis, we will concentrate on the schemes S3 
and S4 of Ref. [7] and the CCS scheme developped in Ref. [8], following the version including quark contributions 
from Ref. [14]. The detailed formulae defining these kernels are explicitely given in the Appendix. We consider 
first the NLO schemes originally defined in Ref. [7] among which only the so-called 53, 54 schemes are valid for 
phenomenological use. We also consider the resummation scheme of Ref. [8], as developed in Ref. [14] by including 



4 In particular, the square root prefactor of the gaussian saddle-point approximation can be merged in the normalization. 

5 The construction of the effective NLO BFKL kernel appears already in Ref. [8] . 
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the quark contribution and denoted CCS. The same procedure can be easily extended to other solutions for NLO 
BFKL kernels proposed in the literature. 

In Fig.l, we show for the three different schemes xnlo{i as a function of 7 for oirg — 0.15 and for different 
values of u. In Fig. 2, one finds the obtained kernels Xeffil, a RG) for the CCS and S3 schemes 6 , after solution of the 
implicit equation (5) for ^(7, ajio). 

In Fig. 3, one displays the parameter-independent values obtained for the set {7c, aRG Xe//(7c) 5 olrg x"//(7c)} as a 
function of aRG- They result from a numerical analysis of the NLO effective vertices (for the CCS and S3 schemes) 
in the vicinity of the saddle-point. Knowing the dependence «rg(Q 2 ), they allow one to predict the behaviour of F 2 , 
up to the determination of the free parameters J\f and Qo, see (8). By comparison we show the corresponding values 
taken by the parameter-dependent set {-y c = i, axLoilc), a XLo(^c)} which are the ingredients of the LO formula (4). 
It is interesting to note that the hard pomeron intercept 0(rg Xeffilc) is compatible with the LO fitted value in a 
physical range of aRG ~ -2 while the LO value a ~ .09 is not physically motivated. 

Let us come now to the quantitative analysis. 

The parameters of the fits are given in Table I. Note that for the LO quantities a is a fitted constant whereas in 
the NLO ones it is given by the standard renormalization group formula (6). However for completeness and inspired 
by theoretical arguments presented in Refs. [10, 11, 13], we have also considered the LO expression (4), with the 
constant a s replaced by the running coupling (6). The corresponding fit parameters are also given in Table I (referred 
as LO'(ai?G))- 

In Fig. 4, we display the results of the BFKL fits to HI data for LO and the S3 and CCS schemes at NLO. As 
is clear from the figure the LO fit is doing a much better job than the considered NLO schemes (even considering 
that the NLO kernels depend on one less parameter). While the qualitative behaviour is correct they fail to take into 
account the quite precise F 2 data, especially at lower Q 2 where they show a too steep behaviour in x B j. Note that 
the S3 scheme is somewhat better as confirmed by the value of the x 2 , displayed in Table I. As obvious from Table 
I, the LO'(«_rg) fit (not represented in the figures) is not successful cither. 

In order to show the behaviour of the different fits with more detail, we display in Fig. 5 the comparison of the 
ratio theory/data for the LO fit (with the expected error bars) and the NLO ones. This figure clearly confirms the 
problems at lower Q 2 . 



IV. ANALYSIS IN MELLIN SPACE 



In this section we want to analyze in more detail the features of the BFKL parametrizations and in particular the 
reasons of the still quantitatively unsatisfactory results of the NLO fits. For this sake, it is important to come back to 
the key ingredient of our analysis, i.e. the dominance of the hard Pomeron singularity expressed by the relation (5). 
As mentionned in section II, this relation is expressed in Mellin space, and our aim is now to make a phenomenological 
test of this relation directly in the ui space, without using the approximation of effective kernels. 

Equality (5) can be checked at NLO using the GRV98 [16], MRS2001 [17], CTEQ6.1 [18] and ALLM [19] parametri- 
sations. These four parametrisations give a fair description of the proton structure functions measured by the HI 
and ZEUS collaborations over a wide range of XBj and Q 2 , as well as fixed target experiment data. The three first 
parametrisations correspond to a DGLAP NLO evolution whereas the ALLM one corresponds to a Regge analysis of 
proton structure function data 7 . Since these parametrisations give a good description of data, we can use them to 
test easily the properties of the NLO BFKL consistency condition (5). This allows us to make a direct computation 
of the Mellin transform F 2 (u, Q 2 ) of the proton structure function and study the NLO BFKL properties in Mellin 
space. Writing the Mellin transform as (cf. formula (1)) 

-/£ (I)' (I0) 

it is easy to realize that 

7WH ™!>, (u) 



6 The 54 scheme gives an effective kernel undistinguishable from S3. 

7 We introduced the ALLM parametrisation in order to avoid biases which could be due to the DGLAP constraints. In fact no significant 
difference appears in the resulting analysis, at least at moderate Q 2 where our analysis is being done. 
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BFKL fit 


a 


Ql 


Af 


X 2 ( (/dof)) 


LO 


0.092 ± 0.010 


0.401 ± 0.012 


.103 ± .002 


1.33 (70) 


LO' (a RG ) 




.055± 0.23 


1.055 ± .032 


8.93 (71) 


NLO (S3) 




3.39 ± 0.23 


.101 ± .003 


3.13 (71) 


NLO (CCS) 




4.27 ± 0.30 


.091 ± .007 


8.62 (71) 



TABLE I: Results of the BFKL fits to the HI data. LO BFKL (a = est.): (lrst line); LO BFKL (a RG (Q 2 )): (2nd line);NLO 
BFKL S3 scheme (3rd line); NLO BFKL CCS scheme (4th line). 



where 7* is the saddle-point for the structure function, or in other well-known terms its effective anomalous dimension . 
Note that under our assumption that the saddle-point of the gluon Green function is transmitted to the full amplitude 
(see previous section) 7* w j c , and thus the relation (11) allows for a phenomenological approach of the consistency 
condition. 

In Fig. 6, we display log F 2 {to, Q 2 ) as a function of log Q 2 for the four parametrisations described above for different 
values of lo. The upper curve corresponds to w = 0.3, and lo varies in steps of 0.1 up to the down curve which 
corresponds to lo = 1.0. The limiting values 8 in to are due to the fact that, one the one hand, we do not want to use 
data at too high XBj to test BFKL NLO properties and on the other hand, the data are lacking at very small xbj- 
We do not see large differences between the parametrisations used in this analysis. The slope in Fig. 6 gives directly 
the value of 7* according to formula (11). The vertical lines define the bins in Q 2 where the slope of F 2 (lo, Q 2 ) can 
be safely evaluated from the curves. 

In Fig. 7, we derive the values of 7* which are the slopes of logi 7 ^ in logQ 2 in the different bins of Q 2 defined 
above, for the four parametrisations. We do not see sizable differences between the parametrisations except at higher 
values of Q 2 , which are anyway not used in our analysis. Hence, in the kinematical range appropriate for BFKL 
phenomenology, the different sets of structure functions, being or not driven by the DGLAP equations, do not give 
rise to noticeable differences. 

After having determined the values of 7*, it is possible to test whether the NLO BFKL formula (5) admits a 
phenomenological verification, i.e. whether -y c w 7*. In other words we have to check the consistency condition 
expressed as 

Xe //(7>, Q 2 ), a RG {Q 2 )) = " . (12) 

ocrgKQ ) 

We considered (12) following the two resummation schemes defined in section II. For scheme S3 x e ^ as a function 
of lo is given in Fig. 8. Note that the different parametrisations agree except at high values of Q 2 , which are anyway 
out of the scope of our BFKL study. 

Interestingly, while the scheme £4 (not shown in Fig. 8) gives the same curves as for 53, the same test could not be 
safely performed for the CCS scheme. The reason is that a spurious pole appears in the quark sector at lo = 1, when 
7* 7^ due to the incomplete momentum sum rule in the quark sector, which gives an effective anomalous dimension 
slightly below the value at lo = 1. This is a technical difficulty which requires a better treatment of the quark and 
gluon sectors simultaneously, which is still a challenge in the NLO BFKL approaches [20] . 

We notice in Fig. 8 that the linear property of relation (12), namely for Xef f(j* i^*, Q 2 ), ctRaiQ 2 )) as a function of 
lo (12) is well verified. We indeed can describe the GRV and MRS parametrisations using a linear fit with a good 
precision. However the predicted zero at the origin lo — is not obtained, even if the value at the origin remains 
small. The fit does not go through the origin and we would need to add a constant term to the linear fit formula. In 
order to quantify the observed discrepancy, we reformulated the phenomenologically obtained curves by the following 
formula: 

X NLO (j*(lo, Q 2 ), a m ) = io/a out , (13) 

where a in = a^Q is the theoretical input (6) while a out is the phenomenologically determined slope of the linear fit 
displayed in Fig. 9. The validity of (12) would obviously require a in = a out . In order to make the comparison, we have 
drawn the straight lines going from the origin, corresponding to the consistency condition (12). 



By various numerical tests, we checked that limiting the Mellin transform range to 0.3 < lo < 1 corresponds to take into account data 
points with 1CT 4 < x B j < 10~ 2 . 
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We show in Fig. 10, the values of on n using the RGE equation (upper curve), and the values of a ou t (lower curve), 
which are always smaller than the values of ajiQ. a out is found to be closer to ai n at low Q 2 but more different from 
this value at higher values of Q 2 . 



V. CONCLUSION 



Summarizing the results of our paper, we have confronted the predictions of BFKL kernels at the level of leading 
and next-leading logarithms (with resummation) with structure function data, using two different proposed types of 
resummation. Our method can be extended to other resummation proposals. 

In a first stage we have proposed to use the "effective kernel" approximation of the NLO-BFKL kernels which, 
associated with the usual saddle-point approximation at high rapidity and large enough Q 2 , allows one to obtain a 
simple two-parameter formula for the structure function F^- The comparison with the similar 3-parameter formula 
commonly used at LO level shows a deterioration of the fits when using two of the known rcsummed NLO schemes 
and a sensitivity to the different types of resummation. 

In order to look for a more model-independent discussion of the discrepancy between precise data and the formulation 
in xsj-space, we perform an analysis of the kernel properties in Mellin-space. For this sake we find an interval in the 
energy-conjugate variable oj where the different Mellin transformed analyses from data give a definite answer within 
reasonable error bars. In Mellin-space, we find that small but sizeable effects give phenomenological deviations from 
the expected theoretical properties of the NLO kernels. 

One possibility is that the saddle-point approximation we introduced is not valid. This simplicity assumption is 
phenomcnologically motivated by its validity already at LO level. Unknown aspects of the prefactors, in particular 
the non-perturbative ones, could play a role in these deviations. 

One way out is to look for higher order effects which could serve as a guide to improved resummation procedures of 
NLO BFKL kernels, hence, it deserves to investigate the phenomenological virtues of other proposed schemes and/or 
trying to make the suitable modifications to the known ones. Also, more precision in the discussion of NLO-BFKL 
predictions will soon be available with the completion of the perturbative NLO impact factors. 

A possibility is to incorporate subasymptotic effects which go beyond the saddle-point approximation and may be 
computed from the NLO formulae [11, 13, 21]. We expect our methods to be valid when incorporating these effects 
in the phenomenological analysis of F2 data and thus further studies in the direction proposed in the present paper 
seem welcome. 
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1. NLO-BFKL Schemes S3 and S4 of Ref.[7] 



Let us define the function 



4xi(7) 



|^Xo(7) 2 - ^Xo(7) - 6C(3) + - 



TT 2 C0S(7T7) 

sin 2 (777) (1 — 27) 



3+ A + M 2 + 37(1-7) \ 

6+ V + Cl) (3 -27XI + 27) J 



-V/'( 7 )-V"(l-7) 



sin(7T7) 



4^(7) 



where 



_ 11CU 2n f K _Q7 tt 2 10n/ 



(A.l) 



(A.2) 



7 



Ca = 3 and n/ = 3, 



and 



0( 7 ) = £(-!)» 



n=0 



ip(n + 1 + 7) - ip(l) , V(« + 2-7) --0(1) 



(n + 7 ) 2 



(n + I-7) 2 



(A.3) 



X o(7) = 2^(1) - V(7) - - 7) • 
We consider the asymptotic expansion of xi(l) m : 

' 67+13^- 



(A.4) 



-K + 



tt 2 +4^'(1) 6C(3)+ 

7 



+ 



395+71 *4-^»-H +OM -^-(44 + K i+ i))^ 



27 



'if 67+13^ ^ 



1 3 



1 + ^'(1) --C(3) 

36 4 v K ' j 7 2^ v ; 



395 + 71^ ^2 



-5- + °W . 



4 36 4 r w / 72 ' ' 108 

and denote di^ the coefficient of -\ in this expansion. In particular we get 

, 67 13n, ir 2 K 0o 11 n f 

-~36~36Cf + i2 + 4" ;dl - 2 ~ '2C~ A ~ 12 ~ ecf ;dl < 3 



(A.5) 



(A.6) 



We write 



where A, B are constants. The coefficient of the term in a in the expansion of \^ is 



Xi 0) (7) = ~(B+ 1/2*0(7)) (V>'(7) + /(I - 7)) - ^Xo(7) , 
whose asymptotic expansion in is given by 

2 

Xi 0) = -(^ + 7^ + 0( 7 2 ))(^ + 2^(1 - 7) + 0( 7 )) - A(± + 0( 7 2 )) - - 4 - + 0(1) • 



(A.8) 



2 7 



2 7 3 <-y 



Denoting d^ ]. the coefficient of \ , one obtains 



Hence, 



Xi (7) - Xi 0) (7) = (di,2 - 4°2)^ + Ki - <4°1)^ + 0(1) • 



Choosing A and B for eliminating spurious singularities in : 



A = -di,i — — :B = -di 2 . 
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(A.9) 



(A.10) 



(A.11) 



we define 



(7) = X (0) (7, w = "X (1) ) = X (0) (7, w) + «(Xi (7) - Xi° } (7)) ■ 



(A.12) 
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The final NLO-BFKL kernel of 53 is finally given by performing the shift 



X (1) (7) = X (1) '(7-^>). 



1 



(A.13) 



For the scheme 54 of Rcf.[7], instead of (A. 7), one starts with 



x (0) (7 ) = X (0) (7)W = ^(o)) = Xo(7) _ 1 _ _L_ + (1 - aA) 



7 1-7 



+ 



7 + \u + aB 1 - 7 + \u) + aB 



One has 



Xi 0) (7) = -4i + T l 



= -A 



7 1—7 
1 1 



2 Xo(j) + B) ^-L + lr A_ 



,7 1-7 

near 0. The singularity coefficients are now 



2 1 +B ) ( 7 2 + (l-7) 2 



+ 0(1) 



Thus, to cancel the divergences of Xi(7) — X 1^(7)1 one i s led to choose 



— 7j ' ^ = — ^1,2 • 



X ( H7) is calculated as for 53. 



(A.14) 



(A.15) 



(A.16) 



(A.17) 



2. NLO-BFKL Scheme CCS [8] as denned in [14] 

Starting from the NLO-BFKL term xi(7) as m [A.l], one first introduces the DGLAP splitting functions and their 
Mcllin transforms, which are 



Giving 



1 0o 



1 1 

+ 



(A.18) 



[V(2 + w)-V(l)] , (A.19) 



Jo 



dzz»P qg {z) = -J- 



Nt f 1 



2N C \ l + 2 + uj 3 + cj 



one defines 



S~1 

A T (u) = A gg (u;) + -^A qg (u), 



with Cf = {N% — 1)/2N C . One can write the pole structure of Xi(7) as 



(A.20) 



(A.21) 



f \ 1 1 
Xi(7) = -77-0 - 



A T (0) , M0)-^_ f /l , 1 



2 7 3 2(1 - 7 ) 3 7 2 (1 - 7) 



7 1-7 



+ finite, 



(A.22) 
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where F, which vanishes when Nf — causing the simple poles to go away, is 
Following the same extraction of spurious poles as previously, we are led to define 

x,W - *,« - + ^-i + F (l + T L-)-W»- A -f^L , (A , 4) 

a function with no poles at all at 7 = and 7=1. Next, one defines 

1 (\ 1 \ A t (lu) A t {uj)-M- 

X i(7,w) = Xi(7) + 7 ~ W 7 + , ..... ) + + n ~ , : ■ (A.25) 

Finally the full kernel is given by 



7 \7 1 — 7 + w/ 7 2 (1 — 7 + w) 2 



x(7^)-Xo(7^)+^^7^ . (A.26) 



which has to be shifted as in (A. 13). 
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FIG. 1: Xnlo as a function 0/7. The curves correspond to a = cte = 0.15 for two values of u). Dark lines: CCS scheme; 
dashed lines: S3 scheme; dotted lines: 5*4 scheme. 
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FIG. 2: lo as a function of"/ for different values of ana- Left: CCS scheme. Right: S3 scheme, see text, org varies between 
0.1 to 0.24 by steps of 0.02. For each value of u in Eq.(5), two values of 7 are found and reported at each side of the minimum. 
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FIG. 3: Determination of the saddle-point variables as functions of org- Top: the saddle-point value y c ; Middle: the effective 
Pomeron intercept org Xeff't Bottom: the effective diffusion variable org Xeff- The LO fixed values correspond to the fit of 
section III. 
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FIG. 4: Results of the BFKL fits to the HI data. 
CCS: Dotted lines). 



LO BFKL kernel (continuous lines); NLO BFKL kernels (S3: Dashed lines, 
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FIG. 5: Ratio theory/data for the BFKL fits to the HI data. LO BFKL: points with error bars; NLO BFKL Scheme 5*3: Dashed 
line; NLO BFKL Scheme CCS: Dotted line. 
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FIG. 6: log F2(lo, Q 2 ) as a function of log Q 2 . Four different parametrisations of the proton structure function have been 
considered: MRS2001 (continuous), GRV98 (dashed), ALLM (dotted), CTEQ6.1 (dotted-dashed) . 
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FIG. 7: Derivation of the anomalous dimension y*(u),Q 2 ) by computing the slope of logF-2 as a function of logQ 2 . 
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FIG. 8: Test of xi^tQ 2 ) f or scheme S3. The results for the four parametrizations are shown for the different bins in Q 2 , (see 
text). 
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FIG. 9: Phenomenological analysis of the effective kernel in Mellin space. Continuous lines: Xeff- Dashed lines: linear 
approximation for as a function of lu. for the different bins of Q 2 and for scheme S3. Continuous straight lines: consistency 
condution (12), see text 



